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ABSTRACT 

The physical mechanisms of the quasar ultraviolet (UV)-optical variability are not well 
understood despite the long history of observations. Recently, Dexter & Agol presented a 
model of quasar UV-optical variability, which assumes large local temperature fluctuations 
in the quasar accretion discs. This inhomogeneous accretion disc model is claimed to de¬ 
scribe not only the single-band variability amplitude, but also microlensing size constraints 
and the quasar composite spectral shape. In this work, we examine the validity of the inhomo¬ 
geneous accretion disc model in the light of quasar UV-optical spectral variability by using 
five-band multi-epoch light curves for nearly 9 000 quasars in the Sloan Digital Sky Survey 
(SDSS) Stripe 82 region. By comparing the values of the intrinsic scatter crj nt of the two-band 
magnitude—magnitude plots for the SDSS quasar light curves and for the simulated light 
curves, we show that Dexter & Agol’s inhomogeneous accretion disc model cannot explain 
the tight inter-band correlation often observed in the SDSS quasar light curves. This result 
leads us to conclude that the local temperature fluctuations in the accretion discs are not the 
main driver of the several years’ UV-optical variability of quasars, and consequently, that the 
assumption that the quasar accretion discs have large localized temperature fluctuations is not 
preferred from the viewpoint of the UV-optical spectral variability. 
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1 INTRODUCTION 


The flux variability of active galactic nuclei (AGN) and 
quasars has been observed for a long time since its dis¬ 
covery. The variability amplitude in the ultraviolet (UV)- 
optical wavelength range reaches almost an order of magni¬ 
tude over several years, which indicates that the mechanism 
causing the variability in an AGN is the mechanism yxm- 


trolling 

the AGN activity (e.g., Ulrich, Maraschi & Urrv 

1997; 

Nandraetal. 1998; Kawaguchi et al] 1998; Giveon et al. 

1999; 

Hawkins 

2003; Vanden Berk et al. 2004; de Vries et al. 

2005; 

Gaskell 

2008!; Wilhite et al, 20081; Liu et al 

2008 

; Bauer et al. 

2009 

Schmidt et al. 2012; Meusinger & Weiss 

2013; 

Kokubo et al. 

2014 

Ruan et al. 2014). 


There are a number of models attempting to account for 
the AGN UV-optical variability. Several models assume that 
the AGN variability is caused by external (non-A GN) f actors, 


collisions (Torricelli-Ciamponi 

et aTf 200 ( 

I) or multiple 

super- 

novae or starbursts near the 

nucleus 

[Terlevich et al. 

1992 !; 


lAretxaga & Terlevichll 19941 : lAretxaga, Cid Fernandes & Terlevic l 


* E-mail: mkokubo@ioa.s.u-tokyo.ac.jp 


119971 : ICid Fernandes, Sodre & Vieira da Silva! l2000h . However, 
these models generally fail to account for the ubiquity of AGN 
variability that has been confirmed by modern large time- 
domain surveys (e.g., IVanden Berk et all 1 20041) . Alternatively, 
several authors have claimed that the AGN variability is due 
to cha nges in the g l obal mass accret i on rat e in AGN accretion 
discs ( Pere^raetal]_ 20061; Li & Caol 1200 8l : ISakata et al.l l201ll : 
IZuo et alJ 120121 : iGu & Lit 120131) . The variable mass accretion 
rate model seems to explain the large variability amplitude 
and the bluer-when-bright er colour v a riabili t y trend often ob- 


Giveon et al.l 1 99 C J~ Hawkins! 20031; Vande n Berk et al 


Cristiani et al. 19 971: UlmdrMaraschi_&Un 2 



1997 


2004 


2012|JSchmidt_etaF 2012): | Gezari et al! 2013 : Kokubo et alJl2014t 


Gallastegui-Aizpun & SaraiedinJ 20141). However, as shown by 


Kokubo et al.l d2014iT (see also ISchmidt et alJ l2012h . the variable 


mass accretion rate model cannot fully account for the strong vari¬ 
ability observed in UV wavelengths. Moreover, the variable mass 
accretion rate model is not preferred when the large difference 
between the AGN UV-optical spectral variability timescale and the 
sound crossing (and the viscous) timescale of the accretion disc, 
which corresponds to the timescale required for global changes 
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2 Mitsuru Kokubo 


of mass accretion rate within the w hole of the accretion disc, is 
taken into consideration (see e.g. , ICourv oisier & Clavell Il99ll : 

lLaMassa et a lJ2014l:ilJttlev & CasellcJl20l4). 

Recently an alternative AGN variability model was pre¬ 
sented bv lDexter & Aeoli J201 lh : the strongly inhomogeneous ac¬ 
cretion disc model, which assumes large local temperature fluc¬ 
tuations in the quasar accretion discs as the cause of flux vari¬ 
ability. Dexter & Agol’s inhomogeneous accretion disc model is 
motivated by recent numerical sim ulatio n studies of thermal or 
magnetorotational instabilities (e.g., iHirose. Krolik & Blaesll2009l; 
Ijiang. Stone & David 1201 3L and references therein). This model 
aims to describe not only the observed amplitude of AGN UV- 
optical single-band variability, but also the unexpectedly large size 


(see e.i 

g.. Poolev et al. 

2007; Dai et al. 2010; Jimenez-Vicente et al. 

20 ii 

Edelson et al. 

2015p. the excess emission in UV wave- 


lengths observed in a co mposit e Hubble Space Telesco pe spe ctrum 
of qu asars dZheng et alJ 1 19971 : [Kawaguchi, Shimura & Mi neshigel 
|200l|), and the stochastic properties o f quasar light curves 
dKellv. Bechtold & Siemiginowskal 20091: iKoztowski et all l20ld: 


^acLeod_sflaL; 2012l;lzu etalJ l2013l : lAndrae. Kim & Bailer-Jonesl 


20131 : Morganson et al. 2014 . and references therein). The Dexter 


& Agol inhomogeneous accretion disc model has received much 
attention, and several authors have claimed that this model can 
reproduce s everal observed properties of quasar UV-optical vari- 
ability (e.g jMeusinger & Weis dl2013l;lRuan et al]|20l4 ; ISun et all 
120141) . For example. iMeusinger & Weiss 1 20131) discussed that the 


observed anti-correlation between the variability amplitude and the 
mass accretion rate of quasars could be explained by the Dexter 
& Agol inhomogeneous accretion disc model. iSun et a 1. d2014l) 
discussed that the Dexter & Agol inhomogeneous accretion disc 
model seemed (at least qualitatively) to account for the timescale- 
dependent colour variability of quasars. 

iRuan et akl d2014l ) constructed a composite difference spec¬ 
trum using two-epoch quasar spec tra (fo r 604 quasars) from the 
Sloan Digita l Sky Survey (SPSS) l lYork et al.ll2000l) following the 
procedure of Wilhite e t all I2005I) . and compared it with the Dexter 
& Agol inhomogeneous accretion disc model. The model’s com¬ 
posite difference spectra are calculated as follows: 5 000 model 
difference spectra between any two successive time steps are pro¬ 
duced, and then the model’s composite difference spec trum is pro¬ 
duced by taking a geometric mean. In iRuan et alJ l2014h . they ob¬ 
tained a reasonable fitting of the model’s geometric mean compos¬ 
ite difference spectra to the observed composite difference spec¬ 
trum, and concluded that quasar UV-optical variability is mainly 
caused by the large localized temperature fluctuations in the quasar 
accretion discs. 

However, it should be noted that because the quasar accre¬ 
tion discs have a large physical size, temperature fluctuations oc¬ 
curring at different radii of the quasar accretion discs, which have 
different mean temperatures, must be causally unconnected. Thus, 
the superposition of independent flares from the localized temper¬ 
ature fluctuations generally means that the inter-band flux—flux 
correlation observed in the two-band light curve of each indi¬ 
vidual quasar becomes weaker and weaker when we take more 
and more separated wavelength band pairs, which contradicts a 
well-known feature of AGN variability: the UV-optical continua 


(e.a.. Choloniewski 198ll 

Krolik et al. 1991; Courvoisier & Clave! 

199 1|; 

Korista et al. 1995; 

Ulrich. Maraschi & Urr\l 19971; Winkler 

19971; 

Cackett, Horne & Winkler 2007; Gaskell 2008; Sakata et al. 

2010; Kokubo et al. 2014(). According to this consideration, the 


successful fitting to the observed difference spectrum by the Dex¬ 
ter & Agol inhomogeneous accretion disc model obtained by 
IRuan et all 120141) is expected to be obtainable only when the com¬ 
parison between the model and observed difference spectra is done 
in a composite sense, because the model’s predicted weak inter¬ 
band correlation can be smeared out when compositing the model 
spectra. As mentioned in iKokubo et alJ d2014f) . at least qualita¬ 
tively, it seems to be difficult to explain the large coherent inter¬ 
band variation within the UV-optical wavelength range (i.e., the 
strong flux—flux correlation for each individual quasar) by local¬ 
ized flares in the quasar accretion discs. 

In this paper, we give quantitative counterarguments against 
Dexter & Agol’s inhomogeneous accretion disc model from the 
viewpoint of quasar UV-optical spectral variability. We calculate 
the intrinsic scatter (see Section [2] for details) of the inter-band 
linear correlation in magnitude—magnitude space predicted by the 
Dexter & Agol inhomogeneous accretion disc model and compare 
it with the observed scatter quantitatively. As a result of the com¬ 
parison, we show that the Dexter & Agol inhomogeneous accretion 
disc model actually cannot explain the tight magnitude—magnitude 
correlation often observed in quasar UV-optical multi-band light 
curves. We conclude that it is suspicious that the large localized 
temperature fluctuations are the main driver of the quasar variabil¬ 
ity, and consequently, that the assumption of a strongly inhomo¬ 
geneous accretion disc is not preferred from the viewpoint of the 
UV-optical spectral variability. 

In Section[2] we introduce the statistical method we use in this 
work to evaluate the intrinsic scatter of the magnitude—magnitude 
plots of the two-band light curves. We describe the database of the 
SDSS Stripe 82 multi-band multi-epoch light curves of quasars, 
and then evaluate the variability amplitude and the intrinsic scatter 
of the magnitude—magnitude plots of these light curves in Sec¬ 
tion [3] We describe the Dexter & Agol inhomogeneous accretion 
disc model in Section 147X1 and then show the details of the model 
light curve calculations in Section l4~2l We compare the data with 
the model predictions in Section [5] Finally, discussion and conclu¬ 
sions are given in Section[6] 


2 METHOD 


The flux—flux plots, and almost equivalently the 
magnitude—magnitude plots, of quasar two-band simultane¬ 
ous light curves are k nown to be well fit t ed by straight lines 
of 


Hagen-Thorn 

1997; Cackett, Horne & Winkld 2007 

; Sakata et al. 

201qT 2011; 

Liraetal. 201 it Schmidt et ahl 201 2 !; 

Kokubo et al. 


AGN light curve quantifies the colour o f the variable component 
(Hagen-Thorn 1997; Kokubo et al. 2014), and the linear regression 
slope of the magnitude—magnitude plot can be an indicator of 
the colour variability of AGNs d Schmidt et alJ 120121) . In this 
work, as described later, we focus on the intrinsic scatter from 
the regression line in magnitude—magnitude space, which can 
be considered as an indicator of the strength of the inter-band 
correlation; the larger the intrinsic scatter is, the weaker the 
two-band magnitude—magnitude correlation is. Since the values 
of the variability amplitude and the intrinsic scatter evaluated 
in magnitude—magnitude space do not depend on the absolute 
value of the flux (i.e., the magnitude unit is dimensionless), it is 
preferable to use the magnitude—magnitude plot because it makes 
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V std [mag] 

Figure 1. Comparison of the values of u-, g-, r-, and z-band variability am¬ 
plitudes for the SDSS Stripe 82 quasars derived by the two methods. \4td 
is defined as V sU \ = y/T, 2 — £ 2 , where E is the sample standard devi¬ 
ation and £ is the root mean square of the photometric error of a single¬ 
band light curve. V is the point estimate value of the variability amplitude 
derived by the Bayesian regression analyses from LINMIX.ERR for the 
magnitude—magnitude plots of the band pairs of u-g , u-r, u-i, u-z for 
V(u), g-r, g-i, g-z for V(g), r-i, r-z for V(r), and i-z for V(i). The val¬ 
ues of V for adjacent band pairs are scaled by 3 for clarity. Black points 
with error bars are the randomly selected measurement points indicating 
the typical values of uncertainty in V, and straight lines indicate V = 14 t d ■ 


the comparison with model predictions easy. Thus, in later sections 
we use the magnitude unit when quantifying the quasar variability. 

The linear regression line of a magnitude—magnitude plot of a 
two-band quasar light curve with measurement errors on both axes 
and with intrinsic scatter can be modelled as [following the notation 


of i Kell vll2007h 


Vi 

Xi 

Vi 

ti 


a + b^i + 6i 
"T €x,i 

Vi + e y,i 

G (cr^) , e x 


G , t y ,i ~ G (aU 


(1) 

( 2 ) 

(3) 

(4) 


where the measurement values of the two-band magnitudes are ex¬ 
pressed as Xi and i/i, in which i labels the measurement epoch; 
i =1, 2, • • •. G(x 2 ) indicates a Gaussian distribution function with 
the variance x 2 and with zero mean, and a tilde (~) means that a 
variable on the left side is drawn from a distribution function on 
the right side, is the independent variable and r]i is the depen¬ 
dent variable, representing the true values of the two-band magni¬ 
tudes. t X: i and are the random measurement errors on Xi and 
yi, whose variances are and <J 2 } i, respectively, ti represents 
the intrinsic scatter, whose variance is assumed to be constant and 
is denoted as ct^,, which corresponds to the square of the standard 
error of the regression. 

There are several methods for evaluating the linear re¬ 
gression line and the intrinsic scatte r fo r data with mea¬ 
surement err ors o n _both axes dGulll Il989l: iD’Agostinil | 
Kell^j2007|^jHo2g 1 j5ovvj&Xam>|j2010L Wall & Jenkins! 


2005 ; 


20121 ; 


Feigelso n & Joge sh Babul 20121 ; Park et al. 20121 . and references 
ther ein). i Kellvl d2 007l) introduced a Bayesian approach to this prob¬ 
lem. [KejJjyOC^’s method assumes that the measurement errors 
and intrinsic scatter are Gaussian, and the probability distribution of 
independent variable (true value of the ith measurement Xi), de¬ 
noted as p(£i|t/>) where ip is the parameter set describing the distri¬ 
bution, is modelled as a weighted mixture of K Gaussian functions. 

1 Kelly! d2007 ) mentioned that it was useful to model using 

this form because it was flexible enough to adapt to a wide variety 
of distributions. In the particular case of the magnitude—magnitude 
plots for the quasar light curves, however, there is observational 
reasoning for adopting the Gaussian function as the probability dis¬ 
tribution of single-band quasar light curves are known to be well 
modelled by a first-order auto regressiv e Gaussian process known as 
a damped random walk (e.g. jMacLeod et al3l2012L and references 
therein). This implies that p(£i \ip) for magnitude—magnitude lin¬ 
ear regression of quasar light curves can be modelled as a Gaus¬ 
sian distribution (K = 1), if the light curve is sampled for a 
sufficiently long duration of time [i.e., longer than the damping 
timescale of the time series, which is known to be_~200 d in 
the rest frame for quasar variabilit y, fe.g.. iButler & Bloomll201 ll ; 
lAndrae. Kim & Bailer-Jonesll2013l) l. Therefore, in the following 
analyses, we assume K = 1, although the choice of the exact value 
of K actually does not significantly affect the results obtained. For 
K = 1, p(£,i\tl>) can be expressed as: 


p(CilV’) = 


C2v 


, i (6 - uY 

r i exp {~2—^ 


(5) 


where r and p, are the parameters of this model. This informa¬ 
tive prior [equation for sample positions enables Bayesian 
data analyses for linear regres sion of data with measurement er¬ 
rors on both axes i lGulllll989l;[Keiivll2007l) . The IDL routine for 
the Bayesian linear regression analysis discussed in I Kelly! d2007h 
(LINMIX_ERR) has been made available in the IDL Astronomy 
Users Librar 5 Q The prior density distributions of the model pa¬ 
rameters and the data likelihood function assumed in LINMIX.ERR 
can be found in lKellvl d2007l) . 


1 http://idlastro.gsfc.nasa.gov/ 
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Table 1. Properties of the individual filter combinations for the SDSS Stripe 82 quasar light curves. 


Band pairs (s-l) 

z 


<*> 

Aobj 

( A/epoch ) 

(Arest)s 

(Arest)z 

A(Arest) 

u-g 

0.0 - 

1.5 

1.025 

3984 

53 

1754 

2314 

560 

u-r 

0.0 - 

1.5 

1.025 

3977 

54 

1754 

3044 

1290 

u-i 

0.0- 

1.5 

1.026 

3964 

54 

1753 

3692 

1939 

u-z 

0.0- 

1.5 

1.010 

3502 

52 

1767 

4443 

2676 

g-r 

0.0 - 

2.0 

1.280 

6715 

58 

2055 

2704 

649 

g-i 

0.0 - 

2.0 

1.280 

6708 

58 

2055 

3281 

1226 

g-z 

0.0- 

2.0 

1.263 

5787 

56 

2071 

3947 

1876 

r-i 

0.0- 

3.5 

1.558 

8783 

60 

2410 

2925 

515 

r-z 

0.0 - 

3.5 

1.517 

7417 

57 

2449 

3548 

1099 

i-z 

0.0 - 

4.5 

1.547 

7510 

57 

2937 

3506 

569 


A band pair is expressed as s-l, where s is the band with the shorter effective wavelength and l indicates the other one. (z), 7V 0 bp and (AT ep0C h) are the mean 
redshift, the number of objects and the mean number of observed epochs for the SDSS Stripe 82 light curve sample of each band pair (see Section0for the 
details of the sample selection). (A res t)s and (A res t)z are the mean rest-frame wavelengths of the bands, defined as the effective wavelengths of the bands 

divided by 1 + (z). A(A re st) is their difference (in units of A). 


The LINMIX_ERR estimator has an advantage over other non- 
Bayesian methods in that it calculates the posterior distribution of 
the parameters for the given data using a Markov Chain Monte 
Carlo (MCMC) method, and hence pr ovid es well-defined and reli¬ 
able parameter uncertainties (see, e.g. JPark et al . 2gTJ). Therefore, 
we can estimate the intrinsic scatter in the magnitude—magnitude 
plot, denoted as <Ji nt (data) for the observation data (Section[3} and 
as < 7 i nt (model) for the model light curves (Section [4.2b . with their 
uncertainties, using LINMIX.ERR. 

In addition to eri nt , we can also derive the single-band variabil¬ 
ity amplitude of the light curve plotted on the x-axis and its uncer¬ 
tainty by way of the model parameter r appearing in equation 0 . 
By definition, r is basically the same as the often-used variability 
amplitude indicator Vm, defined as V^d = \/E 2 — £ 2 , where E is 
the sample standard deviation and £ is the root mean square of the 
photometric error of a single-band light curve (e.g ., Vaugh an et alj 
120031: ISesar et alJ 120071 : lAi et alJEoi(il: IZuo et al]l 20 l 2 h . In this 
work, we denote r as V ( x ): 

V ( x ) = t, ( 6 ) 


where x indicates the photometric band corresponding to the x- 
axis of a magnitude—magnitude plot (e.g., x = u for u-g, u-r, 
u-i and u-z plots), and refer to V(x ) as the variability ampli¬ 
tude. Fig. Q] compares the variability amplitude of g-, r-, i- and 
z-band light curves of SDSS Stripe 82 quasars (see Section [3] for 
details) derived by the two methods described above, i.e., V© and 
V. In Fig. [T] the values of V are in good agreement with those 
of Vstd. This result validates the use of LINMIX_ERR as a reliable 
tool to quantify the quasar variability amplitude, and eventually its 
magnitude—magnitude correlation. 

It should be noted that the broad-band photometric light 
curves contain the flux variability not only of the accretion 
disc continuum emission but also of the broad emission lines 


1983; O’Brien, Goad & Gondhalekar 19951: Baldwin et al. 

199^; 

Korista & Goad 2001; 

Peterson & Horne 2004; Korista & Goad 

2004; Wilhite et al. 200 

3 ; Czernv et al. 2013 

; Kokubo et al.l 20l4 

Chelouche et al. 2014; Hernitschek et al. 

20ll Edelson et al. 


makes the scatter in magnitude—magnitude space Omt(data) larger 
because of the reverberation nature of the emiss ion line variabil¬ 
ity (e.g.. ISakata et al.ll 201 lb ISchmidt et alj 1 20121) . In other words, 
we should consider <Ti n t(data) evaluated by the broad-band quasar 


light curves as an upper bound of the intrinsic scatter of the 
magnitude—magnitude correlation of the pure accretion disc con¬ 
tinuum emission. Therefore, even when the broad emission line 
variability dilutes the coherent variation of the pure accretion disc 
continuum emission to some extent, we can assess the invalidity of 
the Dexter & Agol inhomogeneous accretion disc model by con¬ 
firming the relation <Tm t (model) > < 7 i nt (data). 


3 DATA 


Quantifying the intrinsic scatter in magnitude—magnitude space 
requires simultaneous multi-epoch multi-band photometric light 
curves for quasars. Hence, we use a database of the SDSS Stripe 82 
multi-epoch five-band light curves for spectroscopically confirmed 
quasa rs from the fifth SDSS qua s ar ca talogue I Schneid er et alj 
l20ld) presented bv lMacLeod et al] J 20 I 2 I) . SDSS Stripe 82 region 
is located in the centre of the three stripes in the South Galac¬ 
tic Cap and has been imaged with u, g, r, i and 2 SDSS fil¬ 
ters simultaneously_about 60 times on average from 1998 to 2007 
dAbazaiian et alj 120091) . The average wavelengths for the u-, g- 
, r-, i- and 2 -bands are 3551, 4686, 6165, 7481 and 893lA, re- 
spectivel>@. We cross-match the Stripe 82 light curve catalogue of 
MacLeod et al] 1 2042j) with a catalogue of quasar properties from 
SDSS DR7 IShen et aljl201 % and use the improve d red shift es¬ 
timates dHewett & Wildlboid) listed in IShen et al . ] d 201 l|y s cata¬ 
logue as the redshift of each quasar. 

For completeness, we use all the five-band light curves. The 
a>axis of the two-band magnitude—magnitude plot ( x-y space) is 
always chosen to be the shorter wavelength band of the two. This 
results in a total of 10 band pairs: u-g, u-r, u-i, u-z, g-r, g-i, g-z, 
r-i, r-z and i-z magnitude—magnitude plots. To avoid the Lyman 
a absorption (< 1216 A) entering into the «-, g-, r- and i-bands, we 
focus only on the redshift range below 2 = 1.5, 2.0, 3.5 and 4.5 for 
the regression analyses of band pairs containing the u-, g-, r- and 
i-bands, respectively (as listed in TableQ}. 

It is known that the SDSS Stripe 82 light curves contain some 
outlying photometric points that are several magnitudes fainter (or 
brighter) than adjace nt values, most of w h ich must have a non- 
physical origin (e.g.. ISchmidt et all l20ld : IMacLeod et al.l l20ld : 


2 http://www.sdss.org/dr7/ 
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Figure 2. Intrinsic scatter crj nt versus variability amplitude V for the SDSS Stripe 82 quasars obtained by the regression analyses in magnitude—magnitude 
space for all 10 band pairs (Section 0- Gray dots are the point estimates of the values of crj nt and V. Black points with error bars are the randomly selected 
points of <Ti nt and V to indicate the typical values of uncertainty in the point estimates. 1 a (68.2 per cent, inner dashed line) and 2a (95.4 per cent, outer 
dashed line) contours, which are calculated from a composite a m [-V posterior distribution constructed by summing each individual posterior distribution of 
the sample, are also shown in each of the panel. We show the range of median zb median absolute deviation of V as red vertical bars on the right-hand side of 
the panels, and for comparison, we also show the median values of the variability amplitude (corresponding to X^td) derived by Zuo et al. 2012 for the g-, r- 
and z-band light curves of the SDSS Stripe 82 quasars as blue vertical bars on the left of the red vertical bars in the panels for V ( g ), V ( r ) and V(i). A thin 
dashed straight line in each of the panel indicates V = cr int . The number of objects and the mean redshift (z) for the sample of each band pair are indicated at 
the bottom of each panel. 
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Meusinger. Hinze & de Hooiil 1201 ll : IPalanque-Delabrouille et all 
201ll:lZuoetal.ll2012h . To eliminate these outlying measurements, 
a mechanized light curve was generated by applying a seven-point 
median filter to each of the light curves of the five bands, and then 
the measurements with a residual between the medianized light 
curve and the photometric data larger^than 0. 25 ma g were removed 
(following the procedures of iSchmidt et alJuOlOl) . Then, we ex¬ 
clude quasar light curves that show no variability, i.e., whose sam¬ 
ple varianc e is s maller than the mean measurement error variance 
(e.g.. IZuo et al]|2ol2h . In addition, since we are attempting to fo¬ 
cus on the properties of long-term and multi-epoch variability, we 
exclude from our sample quasar light curves that have less than 20 
photometric epochs. The mean total photometric points (N e p 0 ch) in 
the resulting two-band light curves are listed in Table |T| Then, lin¬ 
ear regression analyses in magnitude—magnitude space for these 
two-band light curves are condu cted using LINMIX_ER R. The po¬ 
tential scale reduction factor R dCelman & Rubinll 1992j) is used in 
LINMIX.ERR to monitor the convergence of the MCMC to the 
posterior, and we do not use data for which the condition R < 1.1 
is not satisfied within 10 000 iterations of MCMC. The final quasar 
sample size iVobj is listed in Table Q] 

In Fig. [2] we plot the results of the regression analyses for the 
SDSS Stripe 82 data in magnitude—magnitude space; i.e., the in¬ 
trinsic scatter a m t versus variability amplitude V for all 10 band 
pairs. Point estimates of the values of (Tint and V for each individ¬ 
ual quasar are defined as the mean and the standard deviation of 
the posterior distributions of these parameters, la (68.2 per cent) 
and 2 a (95.4 per cent) contours of the a, nt -V relation for the SDSS 
Stripe 82 data, which we use as the observational constraint in later 
sections, are calculated from a composite ai nt -V posterior distribu¬ 
tion constructed by summing each individual posterior distribution 
of the sample quasars. 

A red vertical bar on the right-hand side in each of the pan¬ 
els of Fig. [2] shows the range of the median ± median absolute 
deviation of V. And for comparison, we also show the median val¬ 
ues of the variab ility amplitude (corresponding to V^td) derived by 
IZuo et all 120121) for the g-, r- and i-band light curves of the SDSS 
Stripe 82 quasars as blue vertical bars on the left of the red vertical 
bars in the panels for V(g), V ( r ) and V ( i ). In Fig. [2] we can see 
that th e values of the variability amplitude V s n derived bv IZuo et alj 
l2012h are consistent with those of V derived by LINMIX_ERR. 
Again (as already mentioned in Section [2^, this consistency vali¬ 
dates the use of LINMIX.ERR as a reliable tool for quantifying the 
quasar magnitude—magnitude correlation. 

As indicated in Fig. [2] smaller values of <Ti nt indicate tighter 
linear correlation in magnitude—magnitude space, and larger val¬ 
ues of V indicate a larger variability amplitude. The la and 2a 
contours generally satisfy V amt, which means that the two- 
band light curves of quasars are highly correlated. 


4 INHOMOGENEOUS ACCRETION DISC MODEL 
4.1 Dexter & Agol’s formalization 

Steady-state black hole accretion discs residing in luminous AGNs 
or quasars are thoug ht to be well described by t he standard thin ac¬ 
cretio n disc model (Shakura & Sun vaevlll973l : iNovikov & Thornel 
1 19731) . The emitted spectrum from the standard thin accretion 
disc can be written as a sum of blackbody spectra from the 
surface elements of the disc (e.g., iFrank. King & Rainel 1 19921 ; 
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Figure 3. Top: A single realization of the 10-yr model spectral time series, 
assuming ot = 0.40, n = 576 and T* = 100000 K. Middle: g- and i- band 
model light curves obtained from the spectral time series shown in the top 
panel, assuming z = 0, a g = 0.027 mag and rj, = 0.029 mag. Bottom: 
Model magnitude—magnitude plot of the g-i band pair light curve shown in 
the middle panel. The straight line indicates the lineal' regression line drawn 
using LINMIX.ERR. 


Kato, Fukue & Mineshigel l2008l ): 


F v (v obs ) oc 



1 

e hv/k B T ctt (R ) _ l 


dS 


T e ff(R) 


f3GM B nM\ 1/4 (Rm\ 3/A (, 

V ) \r) 1, 



(7) 



1/4 


(8) 


© 2015 RAS, MNRAS 000.mfl4l 























































Constraints on the quasar inhomogeneous disc model 7 


Table 2. Assumed photometric errors in the model five-band light curves (in magnitude unit). 


Group 

i 

2 

3 

4 

5 

6 

7 

8 

9 

< 7 * u 

0.018 

0.021 

0.028 

0.037 

0.051 

0.071 

0.097 

0.14 

0.20 


0.011 

0.012 

0.014 

0.016 

0.021 

0.027 

0.035 

0.048 

0.068 

(7 j' 

0.0099 

0.011 

0.013 

0.015 

0.019 

0.025 

0.034 

0.047 

0.066 

Ci 

0.011 

0.012 

0.014 

0.017 

0.022 

0.029 

0.040 

0.057 

0.081 

°z 

0.016 

0.020 

0.027 

0.037 

0.054 

0.079 

0.11 

0.16 

0.25 

Per cent 

1 

2 

4 

10 

18 

28 

28 

8 

1 


where h, ks, G, and a are Planck constant, Boltzmann constant, 
gravitational constant and the Stefan-Boltzmann constant, respec¬ 
tively. i^obs is the observed-frame frequency, which is related to the 
rest-frame frequency v as v 0 h s = v/(l + z). Mbh and M repre¬ 
sent the black hole mass and the mass accretion rate, respectively. 
dS = RdcfidR is the surface element of the accretion disc, where 
R and <f> are the radial and azimuthal coordinates on the disc plane, 
respectively. The disc inner radius Ri n is taken to be R m = 3Rs 
where Rs is the Schwarzschild radius Rs = 2GMbh/c 2 . 3 Rs 
corresponds to the innermost stable circular orbit of a non-spinning 
Schwarzschild black hole. The dimensionless radial coordinate x 
is defined as * = R/R m , which means the disc inner radius is 

Xin — 1 • 


/ 3GMbhM \ 1/4 

V 8naRf n J 


(9) 


is the disc characteristic temperature, which is the only parameter 
determining the spectral shape of the emitted spectrum L v from 
the standard thin accret ion disc (Perevra et al. 2006; Kokubo et alJ 
12014 iRuan et al J f20 1 4|) . For the SDSS quasars, T* ranges from 
50000 to 200000 K if the Newtonia n valu e of radiative efficiency 
(e = 1/12) is adopted llKokubo et alJl2014h . 

It should be noted that there are several pieces of evidence 
from the studies of hot-dust-poor AGNs, dust reverberation 
mapping and infrared polarimetry. suggesting that emission from 
quasar accretion discs extends well into near-infrared wavelengths, 
which means that an accretion disc does not truncate even at 
a few thousands of Rs (e. 


e.g., Kishimoto, Antonucci & Blaesl 
; Kishimc 


2005]; Tomita et al J 20061; iKishimoto et alJ 120071 


Honig_&Kishimot3 |2010|^_|Hao_etalJ 2010|; iLira et all 201 ll ; 
Koshida et alJ 20141 ; Oknvanskv et alJ 20141 and references 


2008; 


therein). Therefore, we fix the disc outer radius as * 0 ut = 2 16 for 
convenience of the model calculations. Since emissions from the 
outer region of the quasar accretion disc (e.g., x > 10000) are 
negligible in the UV-optical wavelength range, the exact choice of 
the disc outer radius has little effect on the resulting UV-optical 
spectra;_ 

iDexter & Agoll (1 2011 1) and lRuan et al.l J20I4I) calculated inho¬ 
mogeneous accretion disc model spectra by dividing the disc into 
n zones per factor of 2 in radius. The zones are log-spaced in r and 
evenly spaced in <f>. From aniVx iV grid with inner and outer radii 
of Xin = 1 and x ou t = 2 16 , the parameter n is given by 


n = N 


log(2) 

log(2 16 ) 


x N = 


N 2 

~16 


To cover the whole range of the parameter space suggested by 
iDexter & Agoll d201 ll) sufficiently (i.e., n = 10—1 000), in this 
work we consider the cases n = 16, 64, 144, 256, 400, 576, 784, 
1024, 1296, 1600, 1936, 2304, 2704, 3136, 3600 and 4096. 

In Dexter & Agol’s inhomogeneous accretion disc model, 



Magnitude 



Magnitude 


Figure 4. Mean photometric errors of the SDSS Stripe 82 quasar light 
curves as a function of the magnitude (upper panel), and a histogram of the 
i-band magnitudes of the quasar light curves (lower panel). The photometric 
errors and the percentages of the assignment into nine groups assumed in 
the model calculations (listed in Table [2} are shown as points in the upper 
and the lower panels, respectively (see Section FOl for details). 


disc inhomogeneity is given by adding temperature fluctua¬ 
tions to the effective temperature profile [equation ([8}j for each 
zone independently. The effective temperature of each zone T e g 
(x, (j)) is assumed to undergo a dam p ed random walk (e.g. , 
iKell^j^eclfiold&Siemiginowska 2009 1 _ Kozlowski et alJ l201(t 
MacLeod et alj|20ld ; iHemitschek et alJ 20141) with an amplitude 
(jt (in the unit of logT e ff); in the limit n —> oo, the damped ran¬ 
dom walk model becomes time independent with a log-normal dis¬ 
tribution of disc temperatures in each annulus whose variance is 
(I n 10 or) 2 /2 [eq uation (2) of l Dexter & Agol |20j~il) 1. As assumed 
in iDexter & Agoll d201 ll) and lRuan et alJ j2014l) . the characteristic 
decay timescale of the temperature fluctuations Tdamp is fixed to be 
200 d, although the results are insensitive to the choice of Tdamp 
dDexter & Agolll2()Til) . 


4.2 Simulating the SDSS Stripe 82 quasar light curves 

The Dexter & Agol inhomogeneous accretion disc model is char¬ 
acterized by two model parameters, n and or, which represent the 
number of separated zones in an annulus of the disc and the ampli- 
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tude of temperature fluctuations, respectively. In addition to these 
two par ameters, the spectral shape depends on a physical parame - 
ter, T* JPerevra et al.120061 : 1 Kokubo et alJl20l3 : lRuan et alJl2014h . 
The quasar redshift z determines which of the rest-frame wave¬ 
lengths fall on to the broad-band filters. In summary, there are four 
model parameters in o ur model calculations: n, <tt, T* , and red- 
shift z. [Dexter & Agoll J 201 ll ) concluded that, for quasar accretion 
discs, the two parameters n and ar should be in the range of n = 
100—1 000 and <jt = 0.35—0.50 to fit several observational con¬ 
straints, including microlensing accretion disc size measurements, 
quasar UV spectral shape an d the UV-optical variability amplitude 
(see Section ITll. lRuan et alj ( 120141) showed that the composite dif¬ 
ference spectrum of the SDSS quasars can be fitted by the Dexter 

6 Agol inhomogeneous accretio n disc model with the same pa¬ 
rameter ranges as lDexter & Agoll j201 ll) . We decide to explore the 
parameter space of 16 ^ n ^ 4096 and 0.1 ^ or ^ 0.6, which suf- 
fici ently covers the whole range of the parameter space suggested 
bv lDexter & Agoll d201 lh . 

The parameters T* and z are taken to represent the SDSS 
Stripe 82 data. For the S DSS quasars, the parameter T* ranges from 
50 000 K to 200 000 K dKokubo et ali20l4lRuan et alj20l4) . The 
redshift z is taken to be z < 1.5 for u-g, u-r, u-i and u-z band 
pairs; z < 2.0 for g-r, g-i and g-z band pairs; z < 3.5 for r-i and 
r-z band pairs; and z < 4.5 for the i-z band pair (see Section^. 

For a fair comparison between the Dexter & Agol inhomo¬ 
geneous disc model and the SDSS Stripe 82 quasar light curves, 
we have to calculate mock light curves that adequately model the 
actual SDSS Stripe 82 data. The number of measurement points 
and their sampling intervals assumed in the model light curves are 
determined to make them similar to the average SDSS Stripe 82 
light curve; i.e., we assume that we have 10 yr of light curves 
with 1,0, 1, 3, 3, 4, 4, 17, 14 and 9 measurement points in each 
of the 10 yr (56 points in total), and the sampling interval within 
each of the year s is fixed to 4 d (see the middle panel of Fig. [3} 
(e.g. JSako et al]|2014l) . In the actual SDSS Stripe 82 data, the first 

7 y r are the runs o btained as p art of the SDSS-I Legacy survey 
lAdelman-McCarthv et alj(2007l ). and the last 3 yr correspond to 
the SDSS-II Supernovae project runs jFrieman et alj 20085) . 

Moreover, for a fair comparison between the model and 
the data, Gaussian photometric errors are added to the model 
light curves and the same linear regression method (i.e., 
LINMIX.ERR) is applied to the photometric error-added model 
magnitude—magnitude plots. We calculate the 100 realizations of 
the model five-band light curves for each of the parameter sets (n, 
<xt, T*, z), and then divide them into nine groups and assign the 
Gaussian photometric errors to each of the nine groups as tabu¬ 
lated in Table [2] The photometric errors and the percentages of the 
assignment into nine groups are determined to make them compa¬ 
rable to the histogram of the photometric errors of the SDSS Stripe 
82 quasar light curves. Fig. [4] shows the mean photometric errors 
of the SDSS Stripe 82 quasar light curves as a function of the mag¬ 
nitude. First, nine representative values of the photometric errors 
for the i-band (black points in the upper panel of Fig. ID are sam¬ 
pled according to the histogram of the i-band magnitudes of the 
quasar light curves rounded to a whole number (black points in 
the lower panel of Fig. [4}- The percentages of the assignment into 
nine groups are fixed to those of the i-band values (black points 
in the lower panel of Fig. ID, as listed in the last row of Table [D 
Then, the i-band magnitude histogram is shifted in the magnitude 
direction assuming the m ean quasar spectral index of /„ oc u~ 0 ' 5 
dVanden Berk et al.l200lh . and the photometric errors for the it-, g- 
, r- and z-bands are sampled according to the shifted histograms 


(upper panel of Fig. ID- The values of the photometric errors as¬ 
sumed in the model light curves for each of the nine groups are 
listed in Table |D Each of the light curves is assumed to have fixed 
values of photometric errors; i.e., we ignore the changes in photo¬ 
metric errors within each of the light curves. Finally, we apply the 
same exclusion criteria for the simulated light curves as those for 
the SDSS Stripe 82 data described in Section[D 

It should be noted that the assumed values of the sampling 
numbers and the sampling intervals described above do not have 
significant impacts on the regression results, since the detailed in¬ 
formation of the sampling cadence is essentially abandoned in the 
linear regression analyses in magnitude—magnitude space. More¬ 
over, since LINMIX_ERR includes the proper models of the Gaus¬ 
sian photometric errors, the exact values of the assumed photomet¬ 
ric errors also do not significantly affect the regression results. 

The top panel of Fig.[Dshows a single realization of the 10 yr 
model spectral time series, assuming ar = 0.40, n = 576 and T* 
= 100 000 K. As expected, the superposition of independent tem¬ 
perature fluctuations on the accretion disc makes the spectral shape 
of the model spectra highly variable. The middle panel of Fig. [D 
shows the g- and i- band model light curves obtained from the spec¬ 
tral time series shown in the top panel of Fig. [D assuming z = 0, 
a g = 0.027 mag and ai = 0.029 mag, where a g and cr are the as¬ 
sumed photometric errors. The bottom panel of Fig.JDshows the re¬ 
sulting regression line for the g-i band pair magnitude—magnitude 
plot of the two-band light curves of the middle panel. A quantita¬ 
tive discussion on the linear regression results for various values of 
model parameters are given in the next section. 


5 COMPARISON OF THE INHOMOGENEOUS 
ACCRETION DISC MODEL WITH OBSERVATIONAL 
DATA 

Fig- ID shows the intrinsic scatter a mt versus variability amplitude 
V for the Dexter & Agol inhomogeneous accretion disc model for 
the g-i band pair, for various values of ar in 0.1 increments (upper 
panel), and for various values of n = 16, 144, 1024 and 4096 (bot¬ 
tom panel). Since large temperature fluctuations lead to large flux 
variability, the variability amplitude V becomes larger when ar 
becomes larger. On the other hand, since the sum of a large num¬ 
ber of independent fluctuations lessens the stochastic behaviour, the 
variability amplitude becomes small er when n becomes larger [see 
equation (1) oi l Dexter & Agoll ( 1201 lh l. In the upper panel of Fig.[D 
it is apparent that there are fewer points for smaller values of <tt; 
this is because we have excluded the model light curves whose flux 
variability is smaller compared to the assumed photometric errors 
(see Sections[3land l4.2l >. 

In Fig. |D the intrinsic scatter cr;„t is only slightly smaller than 
the variability amplitude V (i.e., <Ji n , < V), which means that the 
Dexter & Agol inhomogeneous accretion disc model predicts a 
weak magnitude—magnitude correlation. Moreover, as we can see 
in Fig. |D the a mt -V relation for various values of ot and n follows 
a linear track from left bottom to top right, which indicates that 
the weakness of the magnitude—magnitude correlation is kept for 
the whole of the model parameter space considered here. Since (as 
shown in Fig.|D the model values of <7i nt and V for various model 
parameters are confined within a certain locus in <Ji nt -V space, and 
we are not interested in the details of the model parameter depen¬ 
dence of them, hereafter we refer to the ai„ t -V relation of the Dexter 

6 Agol inhomogeneous accretion disc model as the region in the 
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Figure 5. Intrinsic scatter rx ulI versus variability amplitude V for the Dexter 
& Agol inhomogeneous accretion disc model for g-i band pairs (dots in 
different colours). Coloured points with error bars are the randomly selected 
points of cr int and V to indicate the typical values of uncertainty in the point 
estimates. Top: rr int -V relation for various values 0.1 ^ cry rC 0.6 (in 0.1 
increments; in different colours). All redshifts, z = 0.0,0.5,1.0, 1.5 and 2.0, 
are plotted simultaneously, and the other model parameters are fixed to n = 
576 and T* = 100 000 K. There are fewer points for smaller values of cry 
because we have excluded the model light curves whose flux variability is 
smaller compared to the assumed photometric errors (see Sectionsl3land l4.2l 
for details). Bottom: The same as the top panel, but for various values n = 
16, 144, 1024 and 4096 (in different colours). The other model parameters 
are fixed to cry = 0.4 and T* = 100 000 K. A thin straight line in each of 
the panels indicates cr int = V. 


&int-V space defined by all the model values of <Ji nt and V for all 
the values of the model parameters, as described below. 

Fig. [6] shows the intrinsic scatter crint versus variability ampli¬ 
tude V for the Dexter & Agol inhomogeneous accretion disc model 
and the SDSS Stripe 82 data for all of the band pairs, expressed as 
contours; the solid red lines are the lcr (68.2 per cent, inner dashed 
line) and the 2a (95.4 per cent, outer dashed line) contours of the 
model CTint-V^ relation for the whole of the parameter space: 0.1^ 
ay ^ 0.6 (in 0.1 increments), 16 < n ^ 4096, T* = 50 000, 100 
000 and 200 000 K, and the specific redshift range for each of the 


band pairs. The solid green lines are the same as the red lines but 
for the restricted parameter space comparable with Dexter & Agol’s 
constraints on ay and n: 0.35 ^ ay ^ 0.50 (in 0.05 increments) 
and 144 ^ n ^ 1024. The black dashed lines are the same as Fig.[2] 
i.e., the la and the 2a contours of the ai„,-V relation of the SDSS 
Stripe 82 data. 

As expected, the inhomogeneous accretion disc model pre¬ 
dicts that a tighter magnitude—magnitude correlation (i.e., smaller 
values of aint) can be achieved if we take band pairs with shorter 
separations of the wavelengths (see AA rcs t listed in TableQ}, namely 
u-g, g-r, r-i and i-z. Actually, if we only focus on the band pairs 
u-g, g-r, r-i and i-z, the model ai n t-U relation seems to be consis¬ 
tent with the observed ai nt -V relation. Flowever, we should note the 
discrepancy in a^t between the model and the data gets clearer for 
band pairs with longer separations of the effective wavelengths; in 
particular, for the u-i, u-z and g-z band pairs, the model-predicted 
regions cover only a part of the data regions, and the data generally 
show smaller values of a mt than those of the model for comparable 
values of V. 

To see the discrepancy more clearly, in Fig. [7] we show the la 
ranges of V and ai nt for the SDSS Stripe 82 data and for the Dexter 
& Agol inhomogeneous accretion disc model as a function of the 
filter combination, obtained by projecting the two-dimensional dis¬ 
tributions of aint and V (shown in Fig. [6J on to the V - and aint-axes, 
respectively. Since in Fig. [6] and the upper panel of Fig. [7] we can 
see that the variability amplitudes V of the SDSS Stripe 82 quasars 
are well reproduced by the inhomogeneous accretion disc model 
with restricted model parameter space comparable with Dexter & 
Agol’s constraints [this verifies the consis tency between model cal¬ 
culations by us and lDexter & Agoll d201 ll) 1. in Fig.[7]we show only 
the results of the model calculations with the restricted model pa¬ 
rameters. As mentioned above, we can clearly see in Fig.[7]that the 
difference for ai„ t between the data and the model predictions be¬ 
comes larger when the wavelength separation becomes larger, and 
the data show smaller values of ai n t than those of the model. For 
the band pairs u-i, u-z and g-z, aint(model) > aint(data) is satis¬ 
fied for about half of the SDSS Stripe 82 quasars, and this dis¬ 
crepancy cannot be reduced by simply adjusting the model param¬ 
eters in the Dexter & Agol inhomogeneous accretion disc model. A 
Kolmogorov—Smirnov test reveals that the difference between the 
data and the model distributions of ai„t for the u-z band pair is sta¬ 
tistically significant, and confirms that the null hypothesis (the data 
are drawn from the model distribution) is rejected even at the level 
of significance a = 0.001. It should be noted that, as mentioned 
in Section [2] ai„t(data) evaluated by the broad-band quasar light 
curves is contaminated by the flux variation of the broad emission 
lines and the Balmer continuum emission, thus the true values of 
aim (i.e., the intrinsic scatter of the magnitude—magnitude correla¬ 
tion of the pure accretion disc continuum emission) must be much 
smaller than those observed. 

These results verify the intuition mentioned in SectionQ] i.e., 
Dexter & Agol’s inhomogeneous accretion disc model cannot ex¬ 
plain the coherent flux variation within the UV-optical wavelength 
range often observed in quasar light curves. 


6 DISCUSSION AND CONCLUSIONS 

In this work, we have examined the validit y of the inhomoge¬ 
neous accretion disc model presented bv lDexter & Agoll d201 ll) in 
the light of the quasar UV-optical spectral variability. As assumed 
in the Dexter & Agol inhomogeneous accretion disc model, be- 


© 2015 RAS, MNRAS 000. HI ITU 











10 


Mitsuru Kokubo 





| o.i 



« o.i 


o.oi 


o.oi o.i 

°int(u.i)[mag] 


0.01 0.1 

a int (g,z)[ m ag] 


0.01 0.1 1 
o int (i,z)[mag] 




Figure 6. Intrinsic scatter crj nt versus variability amplitude V for the Dexter & Agol inhomogeneous accretion disc model and the SDSS Stripe 82 data for all 
of the band pairs. The solid red lines are the lcr (68.2 per cent) and the 2 a (95.4 per cent) contours of the model <r m rV relation for the whole of the parameter 
space: 0.1 ^ ctt ^ 0.6 (in 0.1 increments), 16 ^ n ^ 4096, T* = 50 000, 100 000 and 200 000 K, and the specific redshift range for each of the band pairs. 
The solid green lines are the same as the red lines but for the restricted parameter space comparable with Dexter & Agols constraints on ctt and n: 0.35 ^ <jt 
^ 0.50 (in 0.05 increments) and 144 ^ n ^ 1024. The black dashed lines are the same as Fig. [2] i.e., the lcr and the 2cr contours of the o" mt -V relation of the 
SDSS Stripe 82 data. A thin dashed straight line in each of the panels indicates crj nt = V. 
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Figure 7. 1<t ranges of V (upper panel) and rrj Ilt (lower panel) for the SDSS 
Stripe 82 data (black) and for the Dexter & Agol inhomogeneous accretion 
disc model with the restricted parameter space (green) as a function of the 
filter combination, obtained by projecting the two-dimensional distributions 
of (Ti nt and V (shown in Fig. [6) on to the V- and (T mI -axcs, respectively. 
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cause the quasar accretion discs have a large physical size, the 
several years’ temperature fluctuations occurring at different radii 
of the quasar accretion discs must be causally unconnected. Ob- 
servationally, this property of the inhomogeneous accretion disc 
model ought to result in a weak inter-band magnitude—magnitude 
correlation within the UV-optica l wavelength range (see, e.g., 
IPerevra et al]|2006l : I Gaskelll 1 201 ll) . We have used a large sample 
of long-term (approximately several years), simultaneous five-band 
light curves of quasars in the SDSS Stripe 82 region to evaluate in¬ 
trinsic scatter cri nt (data) in magnitude—magnitude space, and com¬ 
pared cri n t(data) with the Dexter & Agol inhomogeneous accretion 
disc model predictions of the intrinsic scatter (Ti nt (model). We have 
shown that, especially for the two-band light curves for the u-i, u-z 
and g-z band pairs, (Ti nt (model) is generally larger than aint(data) for 
comparable values of variability amplitude V (Fig. [6] and Fig. [7]l, 
which indicates that the Dexter & Agol inhomogeneous accretion 
disc model cannot explain the tight inter-band correlation often ob¬ 
served in the quasar light curves within the UV-optical wavelength 
range. This result indicates that the local temperature fluctuations 
in the accretion discs are not the main driver of the several years’ 
UV-optical variability of quasars, and the assumption of a strongly 
inhomogeneous accretion disc is not preferred from the viewpoint 
of UV-optical spectral variability. 

This conclusion seems to be inconsistent with that obtained by 
Ruan et~ahl l20l4) . in the sense that they concluded that the Dexter 


& Agol inhomogeneous accretion disc model accounted well for 
the spectral variability of SDSS quasars. Flowever, as we have al¬ 
ready mentioned in SectionQ] the difference bet ween our result and 
that of iRuan et alj d20l4 ) is the expected one: iRuan et all j20 14| ) 
compared a geometric mean composite difference spectrum of the 
SDSS quasars with geometric mean composite difference model 
spectra, thus the model-predicted incoherent inter-band flux varia¬ 
tions were smeared out. Because in our analyses the coherence of 
inter-band flux variations for each individual quasar is compared 
with each individual time series realization of the model, we can 
confidently confirm that it is difficult to explain the coherent flux 
variation within the UV-optical wavelength range in each individ¬ 
ual quasar by the Dexter & Agol inhomogeneous accretion disc 
model. 


Moreover, it should be noted that employing unexpectedly 
higher-temperature flares than considered here cannot settle the 
problem. If the blackbody flares always have a very high tem¬ 
perature 0§> 100 000 K), the resulting UV-optical spectral vari¬ 
ability would show a strong inter-band correlation because the 
Rayleigh—Jeans tail (i.e., a power-law spectrum with a„ = 2) 
would dominate the whole of the UV-optical wavelength range. 
However, it is known that the observed colour of the quasar vari¬ 
able components is, generally, significantly redder than that of 
the Rayleigh—Jeans spe ctrum l lPerevra et al .1120061 ; iKokubo et alj 
120141 : IRuan et al] 120141) . which indicates that high-temperature 
blackbody flares also cannot be the main driver of the quasar UV- 
optical variability. 

Our study clarifies the importance of exploring the variabil¬ 
ity models, which can simultaneously explain the general proper¬ 
ties of the several years’ AGN variability, namely, the large vari¬ 
ability amplitude and the tight inter-band correlation within the 
UV-optical wavelength range. One of the possibilities is the repro¬ 
cessing model, which assumes the AGN UV-optical variability is 
cau sed by r eproce ssing of X-ray or extreme UV emission (e.g., 
iKrolik et al.lIl99ll) . Since the fluctuations of this higher-energy 
photon flux, originating from the innermost region of the accre¬ 
tion disc, propagate from inner to outer radii in the light-crossing 
time, the reprocessing model can explain the several years’ tight 
inte r-band correlation within the UV-optical wavelength range (see 
lUttlev & Casellall20l4 for review). 


The reprocessing model predicts that the flux variation oc¬ 
curs at the UV wavelengths prior to the optical wavelengths with 
a time lag of a few h ours (dependi ng on the size of the ac¬ 
cretion disc) (e.g.. ICollier et alj 119991 ; ICackett. Home & Winklej 
120071) . On the other hand, mass accretion fluctuations prop¬ 
agate from outer to inner radii in the accretion disc, which 
means that the flux variations caused by the mass accretion 
fluctuations occur at the optical (longer) wavelengths prior 
to the UV (shorter) wavelengths (e.g., I Arevalo et alj 120081 : 
lUttlev & Casellall20l4) . Within the UV-optical wavelength range, 
the inter-band time delay of the AGN continuum for several 
AGNs are reported to be as expected in the re p rocessing model 


2006; Czerny & Janiuk 

2007; Cackett, Horne & Winkler 


2007 

Goicoechea et al. 2012 

Chelouche 2013; Lohfinketal. 

2014: 


1 et alj 


though results from different observing campa igns are conflict¬ 
ing (see, e.g.. iNandra et al.l 1 19981 : lUttlevI 120061) . Moreover, sev¬ 
eral authors have pointed out that the X-ray luminosity of AGNs 
is ins ufficient to power the l a rge UV-optical variability amplitude 
(e.g., lAntonucci et aP 1 19961: iGaskell 2008|; lEdelson et alll20l4: 
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lUttlev & Casellal 120141) . Obviously, further observational studies 
are needed to test these models. 

Finally, we should note that in this work we only focus on the 
long-term (i.e., several years) spectral variability of quasars. It is 
possible that the long-term and the short-term (several days) vari¬ 
abilities are driven by totally different physical mechanisms; there 
remains a possibility that the short-term AGN variability (with 
small variability amplitude) may be caused b y some kind of lo - 


Czerny et al. 2008|; 

Voevodkin 2011 

; |Ai et al. 2013|; Choi et al. 

2014!: Graham et al. 

20l! Sun et al 

2014; McHardy et al. 20141; 

Edelson et al. 2014, and references therein). To clarify the true na- 


ture of the AGN accretion discs, it is crucial to obtain the multi¬ 
epoch dense-time sampling and multi-wavelength light curve data 
for a large sample of AGN/quasars. Ongoing and future wide and 
deep time-domain surveys, such as the Palomar Transient Factory, 
the Panoramic Survey Telescope & Rapid Response System, the 
Dark Energy Survey, SDSS-IV Time-Domain Spectroscopic Sur¬ 
vey and Large Synoptic Survey Telescope dlvezic et al.ll2008l ). will 
motivate the multi-epoch and multi-wavelength monitoring cam¬ 
paigns for AGN/quasars. These observations will clarify the nature 
of the short- and long-term X-ray-UV-optical-near-infrared correla¬ 
tions, and consequently, of the AGN/quasar accretion disc physics. 
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